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Abstract 

We introduce a new approach to model and analyze Mobility. It is fully based on 
discrete mathematics and yields a class of mobility models, called the Markov Trace Model. 
This model can be seen as the discrete version of the Random Trip Model: including all 
variants of the Random Way-Point Model [13] . 

We derive fundamental properties and explicit analytical formulas for the stationary 
distributions yielded by the Markov Trace Model. Such results can be exploited to compute 
formulas and properties for concrete cases of the Markov Trace Model by just applying 
counting arguments. 

We apply the above general results to the discrete version of the Manhattan Ran- 
dom Way-Point over a square of bounded size. We get formulas for the total stationary 
distribution and for two important conditional ones: the agent spatial and destination 
distributions. 

Our method makes the analysis of complex mobile systems a feasible task. As a further 
evidence of this important fact, we first model a complex vehicular-mobile system over a set 
of crossing streets. Several concrete issues are implemented such as parking zones, traffic 
lights, and variable vehicle speeds. By using a modular version of the Markov Trace Model, 
we get explicit formulas for the stationary distributions yielded by this vehicular-mobile 
model as well. 
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1 Introduction 



A crucial issue in modeling mobility is to find a good balance between the goal of implementing 
important features of concrete scenarios and the possibility to study the model from an analyt- 
ical point of view. Several interesting approaches have been introduced and studied over the 
last years [31 lll[ [TBI I14j . Among them, we focus on those models where agents move indepen- 
dently and according to some random process, i.e., random mobility models. Nice examples of 
such random models are the random-way point and the walker models [31 El fTTl [TIT [16] which 
are, in turn, special cases of a family of mobility models known as random trip model |14j . 

Mobile networks are complex dynamical systems whose analysis is far to be trivial. In 
particular, deriving explicit formulas of the relative stationary probabilistic distributions, such 
as the agent spatial one, requires very complex integral calculus and/or sophisticated tools 
like the Palm Calculus gl El US ES] . 

A first goal of our study is to make the analysis of such dynamical systems more accessible 
to the (Theoretical) Computer Science Community by adopting concepts and methods which 
are typical of this Community. A possible solution for this issue could be that of considering 
agents that walk over (random) paths of a graph. Nice results on random walks over graphs 
are available [21 [11], however, such models are not suitable to consider agent speed variations, 
crossways, parking zones and other concrete aspects of mobile systems. 

We propose a new approach to model and analyse mobility. This approach is based on a 
simple observation over concrete network scenarios: 

It is not so important to record every position of the agents at every instant of 
time and it thus suffices to discretize the space into a set of cells and record the 
current agent cell at discrete time steps. 

We exploit the above observation to get a class of fully-discrete mobility models based on 
agent movement-traces (in short, traces). A trace is the representation of an agent trajectory 
by means of the sequence of visited cells. Similarly to the Random Trip Model, our mobile 
model is defined by fixing the set of feasible traces and the criterium the agent adopts to select 
the next trace after arriving at the end of the current trace. 

We define the (Discrete) Markov Trace Model (in short, MTM) where, at every time step, 
an agent either is (deterministically) following the selected trace or is choosing at random 
(according to a given probability distribution) the next trace over a set of feasible traces, all 
starting from the final cell of the previous trace. 

It is important to observe that the same trajectory run at different speeds yields different 
traces (cell-sequences are in general ordered multi-sets): so, it is possible to model variable 
agent speeds that may depend on the specific area traffic or on other concrete issues. 

A detailed description of the MTM is given in Section [2l We here discuss its major features 
and benefits. 

Any MTM D determines a discrete-time Markov chain A4t> whose generic state is a pair 
(T,i): an agent has chosen trace T and, at that time step, she is in position T(i). 

We first study the stationary distribution(s) of the Markov chain Mx> (in what follows we 
will say shortly: "stationary distribution of P"). We show evidence of the generality of our 
model, derive an explicit form of the stationary distribution(s), and establish existence and 
uniqueness conditions for the stationary distributions of an MTM. 

Our last result for general MTM is a necessary and sufficient condition for uniformness of 
the stationary distribution of T>. 
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The above results for the stationary phase can be applied to get explicit formulas for the 
stationary (agent) spatial distribution and the stationary (agent) destination one. The former 
gives the probability that an agent lies in a given cell, while the latter gives the probability 
that an agent, conditioned to stay in a cell v, has destination cell w, for any choice of v and 
w. 

The knowledge of such distributions is crucial to achieve perfect simulation, to derive 
connectivity properties of Mobile Ad- hoc NETworkS (MANETS) defined over the mobility 
model, and for the study of information spreading over such MANETS [7} [8| [TO] . 

We emphasize that all the obtained explicit formulas can be computed by counting argu- 
ments (it mainly concerns calculating the number of feasible traces passing over or starting 
from a cell). If the agent's behaviour can be described by using a limited set of typical traces 
(this happens in most of MANETS applications), then such formulas can be computed by a 
computer in few minutes. 

Our MTM model can thus serve as a general framework that allows an analytical study of 
concrete mobility scenarios. We provide two examples that show its power and applicability. 

In the first one, we consider the Manhattan Random Way-Point (MRWP) model [10^ [SJ 
H3] . This version of the Random Way-Point model is motivated by scenarios where agents 
travel over an urban zone and try to minimize the number of turns while keeping the chosen 
route as short as possible. We then implement this model as a specific MTM and we derive 
explicit formulas for its stationary distributions. In particular, we provide the spatial and the 
destination distributions for any choice of the cell resolution parameter e > 0. We observe 
that, by taking the limit for e — > 0, our explicit formula of the spatial distribution coincides to 
that computed by using rather complex integral calculus in [10] for the continuous space-time 
MRWP model (in terms of probability density functions). 

Finally, we give, for the first time, the destination distribution of the continuous space-time 
MRWP model as well. Both these formulas have been recently used to derive the first analytical 
bounds on flooding time for the MRWP model |9j. 

Our approach can make the analysis of complex scenarios much simpler: it is just a matter 
of modelling objects and events as ingredients of an MTM. After doing that, you do not need 
to prove new properties or new formulas, you can just apply ours. 

As a second concrete example of this fact, we consider a more complex vehicular- mobility 
scenario: The Down-Town model where a set of horizontal and vertical streets cross each 
other and they alternate with building blocks (see Fig [2]). Agents (i.e. vehicles) move over the 
streets according to Manhattan-like paths and park on the border of the streets (a detailed 
description of the model is given in Section [4,ip , Different agent speeds and red and green 
events of traffic lights can be implemented by considering different traces over the same street 
path. 

Thanks to a modular version of our MTM model, we are also able to analyze this more 
complex scenario. In fact, the main advantage of our approach is that a given scenario can be 
analyzed by simply modelling objects and events as "ingredient" of an MTM, thus obtaining 
the stationary probability distributions directly from our formulas. 

2 The Markov Trace Model 

The mobility model we are introducing is discrete with respect to time and space. The 
positions an agent can occupy during her movement belong to the set of points 1Z and they 
are traced at discrete time steps. The set 1Z might be a subset of for some d = 1, 2, . . ., or 
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it might be some other set. It is only assumed that 1Z is a metric space. 
A movement trace, or simply a trace, is any (finite) sequence T = {uq,u\, . . . ,Uk) of at least 
two points. When we mention points we tacitly assume that they belong to 1Z. The points of 
a trace are not necessarily distinct. The length of a trace T (i.e., the number of points of T) 
is denoted by \T\ and, for each i = 0, 1, . . . |T| — 1, let T(i) denote the i-th point of the trace 
T. A trace T can be interpreted as the recording of the movement of an agent starting from 
some initial time to: for every i = 0, 1, . . . , \T\ — 1, T(i) is the traced position of the agent at 
time to + i ■ t, where r > is the duration of a time step. 

In our model, an agent can move along trajectories that are represented by traces. For any 
trace T, let T star t and T en d denote, respectively, the starting point and the ending point of 
the trace. Let T be any set (possibly infinite) of traces. We say that T is endless if for every 
trace T G T , there is a trace T'eT such that T' staTt = T en d- 

For any point u £ 7Z, T out (u) denotes the subset of traces of T whose starting point is u. Let 

P(T) = {u | T out (u) ^ 0} and S(T) = {(T,i) \ T G T A 1 < i < |T| - 1} 

A Markov Trace Model (MTM) is a pair V = (T, such that: 
i). T is an endless trace set such that |-P(T)| < oo; 

zi). ^ is a Trace Selecting Rule for T (TSR), that is, ^ is a family of probability distributions 
{ipu}ueP(T) such that for each point u G P(7~), ip u is a probability distribution over T out (n). 
Similarly to the random trip model [14] , any MTM T> determines a Markov chain Mx> = 
(5(T), P[T, ^]) whose state space is S(T) and the transition probabilities P[T, ^} are defined 
as follows: for every T G T, 

- Deterministic-move Rule. For every i with 1 ^ i < \T\ — 1, Pr ({T,i) — > (T, i + 1)) = 1; 

- Next-trace Rule. For every T" in T out {T end ), Pr((T, |T| - 1) -> (T", 1)) = ip Tend (T'); 

all the other transition probabilities are 0. It is immediate to verify that for any s G S(T) it 
holds that ^2 r£ s(T) ^ >r (s ^ i") = 1. 

- Stationary Properties. We first introduce some notions that are useful in studying the 
stationary distributions of a MTM. For any u, v G 1Z, T(u, v) denotes the subset of traces of T 
whose starting point is u and whose ending point is v. A crucial Markov chain, determined by 
a MTM T>, is its Kernel: its states are the cells of P, considered as turn points (also known as 
way-points), and the transition probability Pr (u — > v) equals the probability that an agent, 
lying on u, chooses any trace ending in v. Given a MTM T> = (T, Vf), the Kernel of T> is 
the Markov chain Ker(T>) = (P(T), K[T, ^]) where the transition probabilities K[T, ^f] are 
defined as follows: for every u, v G P(T), 

Pv(u^v) = { ^TeT(u^)MT) HT{u,v)^<D 
\ otherwise 

Observe that this definition is sound since, for every u G P(T), it holds that 

veP(T) veP{T)T£T(u,v) TeT out (u) 

where the second equality derives from the fact that T is endless. Along with a MTM T> = 
(T,^) we will use the following notations. For every u G P(T), A*(n) = J2TeT out (u)(.\T\ ~ 
l)V'n(T'). Observe that Ay(u) can be interpreted as the expected length of a trace starting 
from point u. Define also T m (u) = {T G T | T en d = u}. The stationary distributions of V 
and that of Ker{T>) are strongly related as stated in the following 
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Theorem 1 Let V = (T, ^) be any MTM. The following properties hold, a) A map it : 
S(T) — > M is a stationary distribution of V if and only if a stationary distribution a of 
Ker(T>) exists such that 

V(T,i) e S(T) 7r«T,i» = 1 M T start )^ Tatart {T) 

22u£P(T) a ( U ) A *( U ) 

b) A map a : P(T) — >■ M. is a stationary distribution of KeriV) if and only if a stationary 
distribution it of V exists such that 

VueP{T) a(u) = l — — £ 7r (( T ' 1 )) 

E Ter 7r((r,i)) Te ^ t(M) 

Stationary Distributions: Existence and Uniqueness. 

Corollary 2 Let V = (T, be any MTM. The following properties hold, a) V has always a 
stationary distribution, b) V has a unique stationary distribution if and only if Ker(T>) has a 
unique stationary distribution. 

Next proposition shows that the Kernel of an MTM can be any finite Markov chain. 

Proposition 3 Given any Markov chain M = (S, P) with S C 1Z, there exists a MTM V 
such that Ker(V) = M. 

Let V = (T, be a MTM. For any two distinct points u, v € P(T), we say that u is connected 
to v in T> if there exists a sequence of points of P(T) (zq, zi, . . . , Zk) such that zq = u, Zk = v, 
and, for every i = 0, 1, . . . , k — 1, ^TeT(zi z,+i) ^ZiC^) > Informally, this can be interpreted 
as saying that if an agent is in u then, with positive probability, she will reach v. We say 
that V is strongly connected if, for every u,v G -P(T), ii is connected to v. Observe that if P 
is not strongly connected then at least a pair of points u, v G -P(T) exists such that w is not 
connected to v. 

Theorem 4 If V is a strongly connected MTM then T> has a unique stationary distribution. 

Stationary Distributions: Uniformity. Let V = (T, ^) be any MTM. We say that V is 
uniformly selective if Vu G P(T), ip u is a uniform distribution. We say that V is balanced if 
Vu G P(T), \T in (u)\ = |T out O)| Observe that if a MTM V = (T, has a uniform stationary 
distribution then it must be the case that 15(7^1 < oo or, equivalently, \T\ < oo. 

Theorem 5 A MTM V = (T, ^) has a uniform stationary distribution if and only if it is 
both uniformly selective and balanced. 

Stationary Spatial and Destination Distributions. We use the following notations. For 
any trace T G T and for any u G 1Z, define 

#t,„ = |{* G N | 1 < i < \T\ - 1 A Tii) = u}\ and T»-{TeT | # T ,„ > 1} 

- We now derive the function s[y) representing the probability that an agent lies in point v G 1Z 
w.r.t. the stationary distribution tt. This is called Stationary (Agent) Spatial Distribution. 
By definition, for any point u G 1Z, it holds that 

*(«)= E T(M)= E #T, U -7T«T,1)) 

(T,i)eS(T)AT(i)=v TET U 
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If the stationary distribution n is uniform, then s(u) = (1/\S(T)\) ■ X^TeTu ^t,u- 

We say that an MTM is simple if, for any trace T £ T and it £ 1Z, #t,« ^ 1- Then, if the 

MTM is simple and ir is uniform, then it holds 

- Another important distribution is given by function D u (v) representing the probability that 
an agent has destination v under the condition she is in position u. This function will be 
called Stationary (Agent) Destination Distribution. By definition, it holds that 

s(u) s{u) 

We define T u (v) = \T U H T m (v)\ and T u = \T U \ and observe that, if the MTM is simple and it 
is uniform, then 

K(v) = ^ (2) 



3 The Manhattan Random- Way Point 

In this section, we study a mobility model, called Manhattan Random-Way Point, an inter- 
esting variant of the Random- Way Point that has been recently studied in |10j . 
Consider a finite 2-dimensional square of edge length L > 0. A set A of n independent Agents 
move over this square according to the following random rule. Starting from an initial posi- 
tion (xo,?/o)) every agent selects a destination (x,y) uniformly at random in the square (i.e. 
every point of the square has the same probability to be chosen). Then, the agent chooses 
(again uniformly at random) between the two feasible Manhattan (shortest) paths. Once the 
destination and the feasible path are randomly selected, the agents start following the chosen 
route with speed determined by the parameter v. In the sequel, we assume that all agents 
have the same speed v that represents the travelled distance by an agent in the time unit. 
However, a variable agent speed can be easily modelled and analyzed by considering more 
traces for any source- destination pair. Once arrived at the selected destination, every agent 
independently re-applies the process described above again and again. This infinite process 
yields the Manhattan Random Way-Point. 

We here consider a discrete version of the Manhattan Random- Way Point which is in fact 
a Markov Trace Model. Agents will act over a square cell grid of arbitrary-high resolution 
according to a global discrete clock. Every agent, within the next time step, can reach any grid 
point (that can be also considered as a square cell) which is adjacent to its current position. 
We emphasize that, as the grid resolution increases, the time unit decreases. This scalability 
allows to observe the process at arbitrary-small time unit and space resolution while preserving 
the ability to choice any possible agent speed. 

In order to formally define the MTM, we introduce the following support graph G e (V e ,E e ) 
where V e = {(ie,je) : i, j G {0, 1, . . . , N — 1}} and E t = {(u, v) : u, v E V t A d(u, v) = e} 
where, here and in the sequel, iV = \L/e] and d{-, ■) is the Euclidean distance. 
Now, given any point v £ V e , we define a set C(v) of feasible paths from v as follows. For any 
point u, C(v) includes the (at most) two Manhattan paths having exactly one corner point. 
More precisely, let v = (x,y) and u = (x',y') we consider the path having first the horizontal 
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Figure 1: The spatial density func- 
tion is shown by a gradation of gray 
(black corresponds to the maximum 
density and white corresponds to the 
minimum density). The destination 
probability over the cross of agent 
position (L/3,L/4) is shown in gra- 
dation of blue. 



segment from (x,y) to (x',y) and then the vertical segment to (x',y'). The second path is 
symmetrically formed by the vertical segment from (x, y) to (x, y') and the horizontal segment 
to (a/, y'). Observe that if x = x' or y = y', then the two paths coincides. We are now able to 
define the Manhattan Markov Trace Model (in short Manhattan-mtm) (T e ,^ e ), where 
T e = {T | T is the point sequence of a path in C(v) for some v £ V e }, 

and \£ € is the uniform TSR for %. It is easy to verify the Manhattan-mtm enjoys the 
following properties. 

Observation 6 The Manhattan-mtm is balanced, uniformly-selective and strongly-connected. 
So, from Theorems^ and^ the Manhattan-mtm has a unique stationary distribution and 
it is the uniform one. Moreover, since the Manhattan-mtm is simple, the stationary spatial 
and the destination distributions are given by Eq. s U\ and respectively. 

So, we just have to count the size of some subsets of traces (i.e paths in G e (V £ ,E e )). Due to 
lack of space, all calculations are given in App. |Fj The point (ie,je) will be denoted by its 
grid coordinates (i,j). 

The stationary spatial distribution for (T e ,^ e ) is 

_ 3((4N 2 -6N + 2)(i+j)-(4N-2)(i 2 + 3 2 ) + 6N 2 -8N + 3) 
M^J- (N* - N 2 )(4N - 2) [ ) 



We now study the Manhattan Random- Way Point over grids of arbitrarily high resolution, i.e. 
for e — > in order to derive the probability densitiy functions of the stationary distributions. 
We first compute the probability that an agent lies into a square of center {x, y) (where x 
and y are the Euclidean coordinates of a point in V^) and side length 25 w.r.t. the spatial 
distribution; then, we take the limits as 5, e — > 0. We thus get the probability density function 
of the spatial distribution (see Fig. [T|) 

■<x,y)^^(x + y)~^ I {x 2 + y 2 ) (4) 



This is also the formula obtained in |T0] for the classic (real time-space) MRWP model. 
The stationary destination density function can be computed very similarly by applying Eq. 
[5J This is described in Appendix[Fl The probability f^ XQ ^(x,y) that an agent, conditioned 
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to stay in position (xo,yo), has destination (x,y) is 



2L-x — y 

xa+yo 

iL{L(x Q +y a )-(xl+y%)) 

L-x +y 

4L(L(x a +y a )-(xl+y£)) 

L + x -y 
4L(L(x a +ya)-(xl+y£)) 

+ 0O 

+ OO 

+ OO 

+ OO 

+00 



if x < xq and y < yo 

if x > Xo and y > yo 

if x < xo and y > yo 

if x > Xo and y < yo 

if x = x and y = y 

if x- = xo and y < yo (South Case) 

if x < xo and y = yo (West Case 

if x = xo and y > yo (North Case) 

if x > xo and y = yo (East Case) 



(5) 



It is also possible to derive the probability that an agent, visiting point {xq, yo), has destination 
in one of the last four cases (south, west, north, and east) (see Fig. [T]) 



,south 



.north 

( x o,vo) 



yo(L - y ) 



.west 
D ( x o,vo) 



.east 



x (L - x ) 



4L(xo + yo)-4(xg + yg) ' ^*°»°> W<*» 4L(x + y ) - 4(xg + yg) 

We observe that the resulting cross probability, (i.e. the probability an agent has destination 
over the cross centered on its current position), is equal to 1/2 despite the fact that this region 
(i.e. the cross) has area 0. This is crucial for getting an upper bound on flooding time [9j. 



4 Modular Trace Models 

Defining a MTM whose aim is the approximate representation of a concrete mobility scenario 
might be a very demanding task. We thus introduce a technique that makes the definition of 
MTMs easier when the mobility scenario is modular. For example, consider vehicular mobility 
in a city. Any mobility trace can be viewed as formed by the concatenation of trace segments 
each of which is the segment of the trace that lies on a suitable segment of a street (e.g., 
the segment of a street between two crossings). Moreover, given a street segment we can 
consider all the trace segments that lies on it. Then, it is reasonable to think that two alike 
street segments (e.g., two rectilinear segments approximately of the same length), have similar 
collections of trace segments. This leads us to the insight that all the traces can be defined by 
suitably combining the collection of trace segments relative to street segments. It works just 
like combining Lego blocks. 

In the sequel, we use the term trace segment to mean a trace that is a part of longer traces. 
Given a trace (or a trace segment) T, the shadow of T, denoted by St, is the sequence of 
points obtained from T by replacing each maximal run of repetitions of a point a by a single 
occurrence of u. For example, the shadow of (u, u, v, w, w, w, u) is (u, v, w, u) (where u, v, 
w are distinct points). Given any two sequences of points T and T' (be them traces, trace 
segments, or shadows), the combination of T and T', in symbols T ■ T' , is the concatenation 
of the two sequences of points. Moreover, we say that T and T' are disjoint if no point occurs 
in both T and T' . Given any multiset X we denote the cardinality of X by \X\, including 
repeated memberships. 

A bundle B is any non-empty finite multiset of trace segments such that, for every T,T' € B, 
St = St'- The common shadow of all the trace segments in B is called the shadow of B and 
it is denoted by Sb ■ 

Two bundles B and B' are non- overlapping if Sb and Sb' are disjoint. Given two non- 
overlapping bundles B and B' the combination of B and B', in symbols B ■ B' , is the bundle 
consisting of all the trace segments T ■ T' for all the T,T' with T € B and T' € B' . Notice 
that, since B and B' are non-overlapping, it holds that Sb-b 1 = Sb • Sb 1 - Moreover, it holds 
that \B ■ B'\ = \B\ ■ \B'\. 
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A bundle-path is a sequence of bundles P = (B\, B%, . . . , B}.) such that any two consecutive 
bundles of P are non-overlapping. A bundle-path P determines a bundle Bundle(P) = B\ ■ 
-B2 •••-£?£. Observe that |Bundle(P)| = nf =1 \Bi\. Given a bundle-path P, let P s tart and 
P en d denote, respectively, the starting point and the ending point of the traces belonging to 
Bundle(P). 

A route R is a multiset of bundle-paths all having the same starting point R s tart and the same 
ending point R en d (i.e. there exist points R s tart and R en d such that for every bundle-path P 
in R it holds P start = R sta rt and P end = Rend)- 

Informally speaking, a bundle-path is formed by traces having the same shadow; introducing 
such different traces allows to model agents travelling on the same path at different speeds (in 
this way, it is also possible to change speed around cross-ways and modeling other concrete 
events). Moreover, routes are introduced to allow different bundle-paths connecting two points. 
By introducing more copies of the same bundle-path into a route, it is possible to determine 
paths having more agent traffic. As described below, all such issues can be implemented 
without making the system analysis much harder: it still mainly concerns counting traces 
visiting a given bundle. 

A Route System is a pair 9\ = (B, TZ) where: (i) B is a set of bundles, and (ii) TZ is a 
multiset of routes over the bundles of B such that, for every R G TZ, there exists R' G TZ with 

R start = Rend- 

We need some further notations. Let TZ U be the multiset {R € lZ\R s t ar t = u}. #p,t is the 
multiplicity of trace T in Bundle(P). #p,b is the number of occurrences of bundle B in the 
bundle-path P. Moreover # RtB = ^2 p eR#P,B and # B ,u = J2reB ir~T,u, where the sums vary 
over all the elements, including repeated memberships. Let f^B denote the total number of 
occurrences of points in all the trace segments of B, including repeated memberships, that is, 

#B = ^ #B,u 
u in Sb 

A Route System «K = (B,TZ) defines a MTM D[9t] = (T[9t], tf[9t]) where 

(i) TpH] = {T | 3R € TZ BP £ R : T G Bundle(P)} Notice that TpK] is a set not a multiset. 

(ii) for every u G P(TpH]) and for every T G T[9i] out (u), 

^ ]u(T) = W ^ 1 ^ |BundSp)| 

1 1 Re1i u 1 1 PER. 1 v n 

The above probability distribution assigns equal probability to routes starting from u, then it 
assigns equal probability to every bundle path of the same route and, finally, it assigns equal 
probability to every trace occurrence of the same bundle path. 

The "stationary" formulas for general route systems are given in Appendix [G] We here give 
the simpler formulas for balanced route systems. A Route System D\ = (B, TZ) is balanced if, 
for every u G S, it holds that 
\{R G lZ\R start = u}\ = \{Re1Z\R end = u}\. 

We are now able to derive the explicit formulas for the spatial and the destination distribu- 
tions; observe that such formulas can be computed by counting arguments or by computer 
calculations. 

Proposition 7 Let 9\ = (B, TZ) be a balanced Route System such that the associated MTM 
V[*R] = (Tpt], ^pt]) is strongly connected. Then, (i) The stationary spatial distributions of 
V[*R] is, for every u G S, 

s{u) = i y *b± y mth Ah[m] = y ^ y 
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(ii) the stationary destination distributions t) ofT>[?R] are, for every u,v £ 5, 



1 #B,u 



~ «r,AAJ(Rl H iri XI 



see 1 1 R.eTiAR cnd =v 1 1 

Observation 8 In i/ie formulas of the stationary distributions stated in Prop. [7| the fac- 
tor is i/ie on/y one depending on trace segments. When points are homogeneous square 
cells of size length d > 0, the following important interpretation hold. Given a trace seg- 
ment T and a cell u of T, the ratio vt, u = d/(#T.u T ) (where r is the time unit) can be 
interpreted as the (agent) instanteneous speed at cell u in trace T. We can thus define 
the instantaneous slowness as SLOWt, m = 1/vt, u = (#T,ui~)/d. We thus get (#b,u)/\B\ = 
(d/r) ■ (J2tgb slow t,u)/|-B|- Since d/r is a constant, the factor #b,u/\B\ is proportional to 
the average of the instantaneous slowness in u of bundle B: this average will be denoted as 
SLOW b,u- I n order to compute the stationary distributions in Prop. [?] of a given balanced 
Route System, we can thus provide the values of SLOW b,u rather than specifying all the trace 
segments. 



4.1 Application: The DownTown Model 

We now use the Modular Trace Model to describe vehicles that move over a squared city-like 
support. This support consists of a square of (n + 1) x (n + 1) crossing streets (horizontal and 
vertical) and buildings (where n is an even number). Buildings are interdicted zones, while 
veichles move and park on the streets. Streets are in turn formed by parking and transit cells 
and every transit cells of a given street has its own direction (see Fig.s[2]and[3]). Moreover, a 
parking cell has its natural direction given by the direction of its closest transit cells. 
A vehicle (agent) moves from one parking cell (the start) to another parking one (the destina- 
tion) by choosing at random one of the feasible paths. To every feasible path, a set of traces 
is uniquely associated that models the different ways a vehicle may run over that path: this 
will also allow to simulate traffic lights on the cross- ways. 

Every street is an alternating sequence of cross-ways and blocks. We enumerate horizontal 
streets by an increasing even index {0, 2, . . . , n}, starting from the left-top corner. We do the 
same for vertical streets as well. In this way, every cross- way gets a pair of coordinates (i,j). 
We say that a direction is positive over a horizontal street if it goes from left to right and while 
the opposite direction is said to be negative. As for vertical streets, the positive direction is 
the one going from top to bottom and the opposite is said to be negative (see Fig.s[2]and[3|). 
Then, the blocks of a horizontal street i will be indexed from left to right with coordinates 
(i, 1), (i, 3), (i, 5), . . .. Similarly, the block of a vertical street j will be indexed from top to 
bottom with coordinates (3, j), . . .). 



We now formally introduce the DownTown Route System = {B , 1Z }; let's start with 
the bundle set B D . 

[Blocks.] Each (street) block is formed by 4 stripes of m cells each with indexing shown 
in Fig. [3J Two stripes are for transit while the two external ones are for parking use. The 
parking stripe adjacent to the transit stripe with positive direction is said positive parking 
stripe while the other one is said negative parking stripe. 

For every ^ i,j ^ n such that (i odd A j even ) V (i even A j odd ), Block (i,j) has 
the following bundles. Bundle B^(i,j) whose shadow is the stripe having positive direction; 
Bundle B^(i,j) is symmetric to B^(i,j) for the negative direction; For each parking cell 
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Figure 2: The Down- 
Town model with n = 
10. Directions for 
each block are shown 
in gray. In blue are 
shown three possible 
routes. The starting 
cells are in green and 
the ending cells are in 
red. 
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Figure 3: The cross way 
at and its four adja- 

cent blocks are shown, with 
m = 6. The shadows 
of two transit bundles are 
shown blue. The shadows 

of start bundles, and 

— h, are shown in green, 
the other two cases (i.e., 
+- 1-, and H — ) are symmet- 
ric. The shadows of end 
bundles, ++ and -| — , are 
shown in red, the other 
two cases are symmetric. 
The shadow of B c ' is 
shown in violet, the other 
three straight cross bun- 
dles are symmetric. The 
shadow of B^ ,++ (i,j) is 
shown in orange, the other 
seven turn cross bundles 
are symmetric. 
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of index k, there are four start-Bundles B$ k (i,j), Bg k (i,j), Bg k (i,j), and Bg~£(i,j); For 
each parking cell of index k, there are four end-Bundles B~^, + k (i,j), B^~ k (i,j), B~^~ k (i,j), and 
B^ k (i, j). The shadows of the above bundles are shown in Fig. [3l 

As for the trace segments, thanks to Obs. we only describe the average slowness of the 
bundles. For the sake of simplicity, we assume the latter depends only on the cell positions. 
For any transit cell of index k, the average slowness in that cell is slow(/c), where slow() is an 
arbitrary positive function. The slow() function is the same for all transit cells of any of the 
above bundles. Notice that, by choosing a suitable function slow(), we can implement variable 
agent speeds and simulate traffic lights at the cross-ways. For instance, assume there is a 
traffic light on a cross- way having two possible states (red and green): the 2 states alternate 
at regular fixed time (say 1 minute). Then the slowness of cells along every of the four adjacent 
transit blocks should be an increasing function of the distance between the cell (determined by 
index k) and the crossway. The exact form of such increasing function depends on the average 
traffic over that cell (notice that the ratio between the highest and the lowest slowness in the 
same transit block might be order of hundreds). As for the parking cells, we assume that the 
average slowness is equal to a positive constant pause. 

[Cross- Ways.] A cross-way is formed by 12 cells as shown in Fig. [3j We have two types of 
associated bundles. 

For every ^ i, j ^ n such that (i even Aj even ), we have: The 4 straight bundles B^ ,+ (i, j), 
B c'~(^i)> s c + (m)> andB c~ The 8 turn bundles B^' ++ (i,j), B^' + ~(i,j), B^ 
B c Moreover, B^ ++ (i,j), B^ + ~(i,j), B c ' and B^' + (i,j). The relative 

shadows are shown in Fig. [3J Observe that the first sign indicates the sign of the in-direction 
and the other indicates the out-direction. We assume that the average slowness of cross-ways 
cells is set to a positive constant cross. 

Let us now introduce the set of DownTown routes TZ D formed by combining the bundles 
described above. First of all, every route contains exactly one bundle-path. So we can describe 
the bundle-path. For every pair (c, c') of parking cells not belonging to the same block, there 
is (only) one bundle-path that goes from c to c'. The structure of a bundle-path has a start- 
bundle for c, followed by an alternating sequence of block and cross-way bundles, and finally 
an end-bundle for c', i.e., S, C±,B2, C3, . . . , B^-i, Ck, E with k ^ 1. 

Notice that two consecutve bundles of a bundle-path sequence belong to two adjacent cross- 
ways and blocks, and the bundles of a bundle-path are all distinct. Let us describe the 
bundle-paths whose starting cells belong to horizontal blocks. The case of vertical blocks is 
fully symmetric, let (i,j) be the coordinates of the starting horizontal block and let (k,z) be 
the coordinates of the ending block. 

[Vertical End Block.] The bundle-path is easily determined by the following driving 
directions: go straight down horizontal street i toward cross-way (i, z); turn to vertical street 
z till block (k, z). 

[Horiz. End Block.] This case yields in turn two subcases. 

- Case j 7^ z. If k = i then go straight down to horizontal street i till block (i, z); if k 7^ i then 
go to cross- way (i,j + a) (where a = 1 if i < k and a = — 1 otherwise) and turn to vertical 
street j + a towards cross- way (k, j + a), then turn to horizontal street k towards block (k, z). 

- Case j = z. The resulting bundle-path depends on whether the starting cell belongs to 
positive or negative parking stripe. Go to cross- way + a) (where a = 1 if we are in 
the positive case and a = — 1 otherwise) and turn to vertical street j + a towards cross- way 
(k,j + a), then turn to horizontal street k towards block (k, z). 

Typical examples of the above paths are shown in Fig. [2l 
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Figure 4: Asymptotical behaviour of the 
spatial probability distribution of hori- 
zontal positive transit cells (0 < a < 1). 
The role of vertical coordinate i is almost 
negligible since the only transit direction 
of such cells is the positive horizontal one. 
This representation takes no care about 
the slowness of the cells. Clearly, in ver- 
tical positive cells, the roles of i and j 
interchange. 



Observe there are some cross-way bundles on the border that do not belong to any route. For 
the sake of convenience, for every bundle B E B D , we define 

R£1Z D 1 1 

Since every route contains exactly one bundle-path and a bundle occurs at most once in any 
bundle-path, then a-B equals the number of bundle-paths containing bundle B. In Appendix 
IH1 we compute a-B for each kind of bundle in order to get the stationary spatial distribution 
given by Prop. [7J Due to lack of space, below we give such formulas only for a cell in a 
horizontal transit block. As for parking cell and cross-way cells we provide explicit formulas 
in Appendix iHl 

Let A = AbpK^J/m 2 be the normalization constant with Ab[9v^] defined in PropUJ Let u be a 
transit cell of index k (with any k £ {1, . . . , m}) in the positive transit stripe of the horizontal 
transit block (i,j) (i.e. i even and j odd), (with i ^ {0,n}). Then 

, slow(fc) / ., k , , .. 1 . ..\ 
s(u ) = 1 \a{t,J) + —b{j)-\ c(j)\ where 

A \ TO TO J 

n 

a(i, j) = (n - j)(2nj + j + n - i - 1) + (n - 2)j - - + i - 2, 

b(j) = n(n + 1) + - - 2( ? i + and c(j) = (n+ l)(n + j) + n - 3 

An informal representation of the asymptotical behaviour of the above function is given in 
Fig. H 

Acknowledgements. We are very grateful to Paolo Penna for useful comments. 
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A Proof of Theorem [I] 

In order to prove Theorem [TJ we need some preliminary lemmas. 

Lemma 9 Let T> = (T, V I / ) be any MTM and let n be a stationary distribution of T>. Then, the 
following properties hold. 

(a) For every (T, i) G S(T), 

n((T,i)) = n((T,l)) 

(b) For every T G T , 

7r«T,l» = ^w(T) 7r (( T '- 1 )) 

T'eT^(T start ) 

(c) For every u G P{T), 

Y <(t,i)) = y ^« T ' 1 )) 

TGT out («) TeT in (n) 

Proof. Since 7r is a stationary distribution of T>, it holds that, for every (T,i) G ^(T), 
n({T,i)) = Y K((T'j))Pr({T',j)^(T,i}) 

(T',j)£S(T) 

If i > 1 then Pr {(T',j) ->■ (T, i)) = for any T' ^ T or j ^ i - 1. Thus, for every 1 < i < |T| - 1, 

7r«T,i}) - 7r((r,t-l» 
This implies that, for every (T, i) G S(T), 

7r«T,<»=7r((T,l» 

and property (a) is proved. 

If instead i = 1 then Pr((T',j) (T, 1)) = whenever T^ d ^ T stort or j ^ |T'| - 1. Thus, 

7r«T,l» = Y *{{T',\T'\-m T3taTt (T) 

T> £T"HT start ) 

= i>T atart (T) Y <( T '^)) (from property (a)) 

T'£T™(T start ) 

This proves property (b). 

From property (b), it derives that 

Y <(t,i)) = Y E A{T',\r\-m u {T) 

TeT out O) TeT° ut (u)T>eT in (u) 

= E ^{T) E <{T',\T'\-l)) 

TeT out («) T'GT in (u) 

= 7r((T',|T'| - 1)) 

T'er ta (u) 

53 7T((TM» 

and property (c) is proved. □ 
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Lemma 10 Let V = (T, v f') be any MTM and let a be a stationary distribution of Ker(D). Let 
tt : 5(7") — > R be defined as follows: 

V(T, i) e 5(70 7r«T, %)) = = 1 M T start )ifr. tart (T) 

E« e p(7-) CT ( U ) A *W 

TTien, n is a stationary distribution ofV. 

Proof. First of all, we verify that 7r is a probability distribution over 5(7"): 

£ 7r«T,»» = ]T = Lu«( U ^ (Trtart) ^'« rt(T) 

E ^(T start )VT start (T) 

(Ta)eS(T) 

\T\-\ 

E E E *(«)^cn 

tieP(r)TGT out (ti) »=i 

2 ^(«) E (m-i)^(T) 

-iGP(T) TeT° ut («) 

ct(u)a*(m) 

^eP(T) 



For proving that 7r is a stationary distribution of V it suffices to show that, for every s G 5(7"), 

7r(s) = ^ 7r(r)Pr (r — >• s) 

reS(T) 

For every s S 5(7"), let A(s) = Eres(T) 7r ( r )P r (r — >• s). Let (T, i) be any state in 5(7"), we distinguish 
two cases. 

(1 < i < |T| - 1): In this case Pr ((T',j) -> (T, i» = for any T' ^ T and Pr ((T, j) -> (T, i» = for 
any j ^ i — 1. Thus, 

T'|-l 

A (M) = E E A(T',j))Pr({T',j)^(T,i)) 

T'ET J = l 
\T\-1 

= K((T,j))Pr((T,j)^{T,i)) 

3 = 1 

= 7r((T,*-l»Pr((T,*-l>->(T,i)) 
= 7T«T,i-l» = 7r((T, i)) 

where the last equality derives from the definition of tt. 
(i = 1): In this case Pr {(T',j) ->■ (T, 1)) = for any V and j such that T' end ^ T start or j ^ |T'| - 1. 



E u gP(T) 


ct(u)A*(u) 




1 


J2uGP(T) 






1 


EmgP(T) 


<r(u)A*(tt) 




1 


E u gP(T) 




1 





2 



Thus, 



|T'|-1 

A«T,1» = E E *({T',j))Pr((T',j)^{T,l)) 

T'eT j=l 

E n((T',\T'\-l))Pv((T',\T'\- 1) <T, 1» 

T'er in (T start ) 

E ^(TMT'l-l))^^^) 

T'GT in (T start ) 
= ^T start (T) E 7T«TMT'|-1» 

T'er»>(r« tart ) 
T'eT in (T star t) e ('> 

E E -(«)<•••« en 



E Me P(T)<^)A*(«) u£P(r)T , er(UiTstart) 

Observe that Et'€T(« T start ) ^u{T') is equals to Pr (u — > T start ) (i.e., the transition probability 
from state u to state T star t of the Markov chain Ker(T>)). It follows that 

A«T,1» = -, x E ^(«)Pr(«-^T rtort ) 

Since er is a stationary probability distribution of Ker(D), it holds that 

E er(w)Pr ( U T start ) = v{T s tart) 
u£P(T) 

Hence, 

A((T,1)) = - r a(T stort ) = 7r«T,l» 



□ 



Lemma 11 Let 2? = (T, &e any MTM and let it be a stationary distribution ofT>. Let a : P{T) 
be defined as follows: 

VueP(T) a(u) = l — — E 7r (( T ' 1 )) 

Then, a is a stationary distribution of Ker(T>). 

Proof. First of all we prove that a is a probability distribution over P(T): 

E '(«» - E Er€ ^ (( T,i» E * (<IU » 

ueP(T) ueP(T) ^ T£T u ' /; rer™'(«) 

= y i (r d) E E 7r ^ T ' 1 )) 

= 1 
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For proving that a is a stationary distribution of Ker(T>) it suffices to show that, for every u S P{T), 

a(u) = cr(i!)Pr — > u) 

veP(T) 

We use the following abbreviations <;(u) — J2 v eP(T) CT ( w )P r ( v ~> u ) an d a = y] \((T l)) ' ^ holds 

that 

E a E tt((T, l»Pr(v-)-u) 

v€P(T) TeT oM (v) 

= « E E 7r((T,l»Pr 

•uSP(T) T6T™'(u) 

= « E Pr H«) E 7r « T > 1 )) 
= « E ( E ^cn) E 7r « T ' 1 )) 

v£P(T) \TeT(v,u) J TeT out (v) 

= a E ( E ^Cn) E 7r (( T '' 1 )) (by Lemma [H(c)) 

= « E E E T«rM))^(r) 

ueP(r) TeT(v,u) T'er in (») 

= a E E 7r (( T ' 1 )) (by Lemma [5] (b)) 

ueP(r) TeT(v,u) 

= a E ^« T ' 1 )) 
rer in (u) 

= a (by Lemma [3(c)) 

= cr(«) 

□ 

Proof of Theorem [H Consider property (a). The "if" part is equivalent to LemmafTUI Now we prove 
the "only if" part. Let ir be a stationary distribution of T>. From Lemma [TT1 the map a : P(T) — > K 
defined as 

V.eP(T) .( U ) = J E 
is a stationary distribution of Ker(T>). For every (T, i) £ S(T), let 

L lie p(r)' J ( tl ) A *( u ) 

From this definition and from Lemma [S] (a), for showing that A = tt it suffices to prove that, for every 
T e T, A((T, 1)) = tt((T, 1)). It holds that 

A (( T ,!)) = ^ w 1 M T start)^T start ( T ) 

L, ue p(T) cr ( u ) A *( u ) 

-^(T s tart)lpT etart (T) 



J2u£P(T) £ r , er i((r',i» ET'er° u *(«) 7r (( T '' 1 )) A *( M ) 

_7rf7T'.lVl 

-a(T start )ip Tgtart (T) (6) 



E T , er7 r((T',l)) 



E„eP(T) Er'e7*»*(u) ^(C 7 "- !)) A *( w ) 
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Now, observe that 

E E T«TM»A*(to = E E ^M)) E cm -iM.cn 

u£P(T)T'eT° ut {u) ueP(T)T'eT mt (u) T"eT out («) 

= E E (|r"|-i)^(T") E 7r (( T '' x )) 

uGP(T) T"6T° ut (u) T'eT° ut (u) 

Since 7r is a stationary distribution of I?, from Lemma[5](b) and (c), it holds that (since T" E T out (u)) 

MT") E ^((TM)) = 7r((T",l» 

T'6T° ut («) 

Thus, 

E E 4(tm))a*( U ) = e E (m-iM.cn e 

uSP(T) T'eT out (tt) u£P(T)T"eT<" lt (u) T'eT out (ti) 

= E E (|T"|-l)7r«T",l)) 

uGP(T) T"6r™'(«) 

|T"|-1 

= E E E ^(( T "' i )) (by Lemma I (a)) 

uEP(T) T"GT out (u) i=l 

= e ^« r ">*» 

<T",i)er 

= 1 

By combining this with Eq. [6l we have that 

A«T,1» = E <( T '^)H T start)^T atart (T) 

T'eT 

= ^T start {T)a(T st art) E <( T '^)) 



T'eT 

= i>T Btart (T)- — 1— — e 7r « T ". 1 )) E ?r (( T '' 1 )) 

= fe tart (T) E 7r (( T "' !)) 

T"6T™'(T, tart ) 

= tt((T,1)) (by Lemma El (b) and (c)) 

Now consider property (b) . The "if" part is equivalent to Lemma 1111 It remains to prove the 
"only if" part. Let a be a stationary distribution of Ker(V). From Lemma [TOl the map it : S(T) — >• R 
defined as 

V{T,i)€S(T) n((T,i)) = 1 M T 8tart )^ Tatart {T) 

is a stationary distribution of V. For every u £ P(7~), let 



- E, gr V,D) S ^1» 
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We will prove that ^ — a. First of all, observe that 

£>((T,1» = ^ ( \K < M T start)^T start (T) 



Thus, we have 
?(«) 



ueP(T) 
1 



XLeP(T) 


ct(m)A^(u) 




1 


SueP(T) 


a(u)A-qr(u) 




1 


XLeP(T) 


<t(m)A$(u) 




1 


Z^eP(T) 


(t(m)A^(u) 




1 


S«eP(T) 





E °( T start)^T start { T ) 
TET 

E E '("^cn 

tiGP(T) TeT° ut (") 

E E 

E (since 

E MT) = 1) 
(since E] 17 ( M ) = 1) 

uGP(T) 



Erer 4(T,i)) E 

E <r(u)A*(u) E 7r (( T ' 1 » 
msp(T) rer°°'(") 

E 0-(u)A*(u) E ? M T start)lpT start (T) 

ueTtr) tet-M E ^ p ( r ) a ( u ) A *^) 

E °( T starMT Btar A T ) 

TeT out («) 

E CT W«( T ) 

T67 -out( u ) 

E ^u{T) 
TeT° ut (u) 

a{u) (since E W T ) = X ) 



□ 



B Proof of Corollary [2] 

Since Ker(T>) is a finite Markov chain, property (a) is a consequence of Theorem [T] (a) and a standard 
result that says any finite Markov chain has a stationary distribution. Although it is a standard result, 
it is not so easy to find a clear reference for it. Thus, we outline a proof of this result. Let M. — (S, P) 
be any Markov chain with 151 < oo. Let n — \S\ and let P = (pij)i,jeS- Consider the subset of the 
n-dimensional Euclidean space U = {(x±, X2 ■ • ■ , x n ) Vi xi ^ A Y^i=i x i — !}■ L : U — > U be the 
map defined as follows: 

n n n 

Yx = (x 1 ,...,x n ) l(x) = cy^Pjix^y^p^xj, . . . ,y2pj n xj) 



3=1 



3=1 



' / j 1 
3=1 



It is easy to verify that, Vx G U, L(x) <E U. Indeed, 53j=i Pji x j ^ for all i. Moreover Y^i=i Sj=i Pji x j ' 
y~lj—i x j Ylj— i Pji — Yl^—i x j — 1; since X)"=i Pji = 1 f° r au Observe that any fixed point of L, that 
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is, a point u G U such that L(u) = u, is a stationary distribution for the Markov chain A4. The 
existence of at least a fixed point is guaranteed by Brouwer's fixed point theorem (see for instance 

my- 

Every continuous map from a convex compact subset if of a Euclidean space to K itself 
has a fixed point. 

Indeed, it is easy to see that U is compact and convex and L is a continuous map. 
Property (b) immediately derives from Theorem [1] (a). 

□ 



C Proof of Proposition [3] 

Let P = (puv)- Define T — {(u,v) | p uv > 0}. Observe that for every (u,v) £ T, there exists at 
least a w G S such that (v, w) G T, since ~Yl x ^sPvx = This implies that T is endless. Also notice 
that P(T) — S. For each u S P(T), define ip u (u,v) — p uv for every (w, u) G T out (w). Clearly tjj u is 
a probability distribution over T out (u). Thus, ^ = {Vv}mgp(T) i s a TSR for T ■ As a consequence 
V = (T, *) is a MTM. Since for every u, v G P(T), EtgT(« «) ^«( T ) = 11 holds that K er{V) = M. 

□ 



D Proof of Theorem [4] 

Lemma 12 7/ 2? is a strongly connected MTM then Ker(T>) is an irreducible Markov chain. 

Proof. Recall that a Markov chain is said to be irreducible if it is possible to go from every state to 
every state (not necessarily in one step). Consider any two states u, v G P(T) of Ker(T>). Since V is a 
strongly connected, there exists a sequence of points of P(T) (zq, Z\, . . . , Zk) such that zq = u, z^ = v, 
and, for every « = 0, 1, . . . , fc — 1, 

Pr{z t ^z l+1 ) = Yl ^( T ) > 

T£T{zi,Zi + l) 

Thus, 

fc-1 

Pr(u^v) {k) ^ Yl Pr (z t ^ z t+1 ) > 

t=0 

where Pr (u —> v) denotes the probability of going from u to v in k steps. Thus, there is a positive 
probability of reaching v from u. 

□ 

Proof of Theorem [4j A standard result says that a finite irreducible Markov chain has a unique 
stationary distribution (see for instance [T5] Chapter 11 - Theorem 11.10, here irreducible chains are 
called ergodic chains). From this result and Lemma[T21 it follows that Ker(T>) has a unique stationary 
distribution. Then, from Corollary [5] (b) , also T> has a unique stationary distribution. □ 



E Proof of Theorem [5] 

We first need the following 

Lemma 13 Let T> = (7", \&) be any MTM that is both uniformly selective and balanced, then Ker(T>) 
has the following stationary distribution a: \fu G P(T)> cr(u) = i-y-j"^ 
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Proof. It holds that, for every u € P(T), 

veP{T) veP(T) TeT(v,u) 



\T(v u)\ 

}^ o-(v)— — ' (since V is uniformly selective) 

v |r out (^)l \T(v,u)\ 
2^ m |r° ut (w)i 



E 



|T(«,«)| 



171 

ir in (u)| 

- (since T> is balanced) 



171 

|T out (u) 



iri 



Proof of Theorem [5j Let 7r be a uniform stationary distribution of V. Then 

1 



□ 



V(T,i) 7r((T,<» 



\S(T)\ 



From Lemma [5](c), it immediately derives that T> is balanced. The uniformly selectiveity derives from 
Lemma [5] (b). 

Conversely, assume that T> is both uniformly selective and balanced. From Lemma fT3l the following 
probability distribution a over P{T): 

Y7 c P m < \ \r out (u)\ 



\r\ 

is stationary for Ker(V). Thus, by Lemma HU1 the following map 7r is a stationary distribution of T>: 

V(T,i) 6 S(T) n((T,i}) = = 1 ^ (JW t )fc ta , f (r) 

Firstly, observe that the following holds: 

£ (t(u)a*(u) = X! l -^wr (iri-i)^u(T) 

«sp(T) ueP(T) ' ' TsT° ut («) 

= y rC!M y (|T| - 1) 1 

Zv 7- w 1 ; |r out (^)l 

= E E (m-D 

|7"| 
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Then, it holds that, for every (T,i) G S(T), 



n((T,i)) = ^ —. r-r , M T start )tl) T3tart {T) 

y ' r tT(T start )^T. tort (r) 



|5(T)| 

m |r out (T start )| i 

|S(T)| |T| |r out (^*art)| 

i 



l^(T)| 

Thus, 7r is uniform. □ 



F The Manhattan-mtm: Stationary Distributions 

Let's remind the key-property that allows to compute such stationary distributions. 

Observation 14 The MANHATTAN-MTM is balanced, uniformly- selective and strongly- connected. So, 
from Theorems^ and^ the Manhattan-mtm has a unique stationary distribution and it is the uni- 
form one. Moreover, since the MANHATTAN-MTM is simple, the spatial and the destinational stationary 
distributions are given by Eq. s [7] and fJl respectively. 

Since the length of a path from a point to a point {if , j') is \i — i'\ + \ j with some calculations, 
we obtain 



\S(T e )\ 



(7) 



( N 4 _ N 2 )(AN - 2) 

Now let rj(u,v) be the number of paths starting from u and visiting v. As for the Manhattan-MTM 
we get 



rie{ii',j')S,3)) 



2N - 


-i-j 


if i' 


< i 


and j' 


< j 


i+j 


+ 2 


if i' 


> i 


and j' 


>3 


N + 


1 - i + j 


if i' 


< i 


and j' 


>J 


N + 


1 + j 


if i' 


> i 


and j' 


<J 


N 2 - 


-Nj 


if i' 




and j' 


<J 


N 2 - 


-Ni 


if i' 


< i 


and j' 


= 3 


N + 


Nj 


if i' 




and j' 


>J 


N + 


Ni 


if i' 


> i 


and j' 


= 3 


2N 2 


~2N + 1 


if i' 


= i 


and j' 


= 3 



Thus 



r ^3) = E E Ve((i'J%(i,j)) 

0<i'<N-l 0<j'<JV-l 

= (4N 2 — 6N + 2)(i + j) — (4N - 2)(i 2 + f) + 6iV 2 



8A^ + 3 



From Eq.s[TJ |9l [3 we get the spatial distribution 

3((4iV 2 - 67V + 2){i + j) - (4N - 2){i 2 + j 2 ) + 6N 2 -8N + 3) 



(iV 4 -iV 2 )(4iV-2) 



(8) 



(9) 



(10) 



Our next goal is to study the Manhattan Random- Way Point over grids of arbitrarily high resolution, 
i.e. for e — ¥ 0. So, we will need to derive the probability- densitiy functions of the stationary distribution. 



9 



Let us compute the probability that an agent lies into a square of center (x, y) (where x and y are the 
Euclidean coordinates of a point in V e ) and side length 26 w.r.t. the spatial distribution. 

Then the probability density function of the spatial distribution is given by 



f(x, y) = lim lim — fs >e (x, y) 



jp(x + y)- ^ i (x 2 +y 2 ) 



(11) 



We now compute the agent- destination distribution. This implies that the number of feasible paths 
visiting point u and ending in point v equals the number of feasible paths starting from u and visiting 



v, i.e. 



T e,(i ,j )(iJ) = Ve((i,j), (iodo)) 
We now replace Eq .s [Inland 151 into Eq. [2], and get the destination distribution 

Ve((i,j), (io, jo)) 



o)(*»j) 



(12) 



(13) 



r e (io,io) 

Let us now compute, the probability that an agent, visiting point (xo,yo), has destination lying into 
the square of center (x, y) and side length 26 (where (^OjJ/o) arL d (x,y) are the Euclidean coordinates 
of points in V e ). By definition of d e , u (v), it follows that 



f( Xo ,y ),5,e( X >y) = E 

«e{f(x-5),4 



E 



r(y+5)} 

The probability density function of the destination distribution is 

1 



e,(£*o,$ito) 



f(x , yo )(x,y) = \im\im-^f {s:0tyo)Ae {x,y) 
By combining Eq.s[l3l 13 [9j and some calculations, we get 



f(x , vo ){x,y) = < 



2L-xq-yg 

iL{L(x a +y a )-(xl+yl)) 

xp+yo 

±L(L(x +y a )-(xl+yl)) 

L-xg+yg 

iL(L(x +y a )-(xl+yl)) 

L+xg-yg 

iL(L(xg+y a )-{x +yl)) 

+oo 
+oo 
+oo 
+oo 
+00 



if x < xq and y < yo 

if x > xq and y > yo 

ii x < xq and y > yo 

if x > xo and y < yo 

if a; = xo and y < y (SOUTH Case) 

if x < x and y = yo (WEST Case 

if x = x and y > y (NORTH Case) 

if x > xo and y = yo (EAST Case) 

if x = xo and y — yo 



(14) 



We now want to compute the stationary agent-destination distribution where the probability den- 
sity function does not exist. To this aim, for every segment s G { SOUTH, WEST, NORTH,EAST }, 
we evaluate the probability that an agent, visting point (xo,yo)i bas destination lying into the subseg- 
ment of s, having any given center and length 26, conditioned to the event that the destination belongs 
to segment s. From Eq. [JJ we observe that, for any segment s, it holds that 



Q u {v\v e s) = 



r u (v) 



We apply the above equation to each of the four cases and get, respectively, 



south ( , _ ^{(n/l)( v - S ),(n/l) {v+S ) } vMN/L)xq,j), ((N/L)xp, (N/L)y )) 
/(%9oWfW E m o,(N/L )vo} vMN/L)xo,j),(WL)xo,(N/L)yo)) V J ° 
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/•West l \ _ ^^xy^i^)y^-u,,yi-,,±,,y^-ru lt ■ , ■ , , - ■ . . f ^ 

^°<y°^ [ > J2 ie{0 , {N/L)xo}Ve ((i,(N/L)yo),((N/L)x ,(N/L)y )) 



T. l e{(N/L)(x-S)XN/L)(x+s)} Ve((i, (N/L)y ), ((N/L)x , (N/L)y )) 



f north 



T, 1 e{(N/L)(y^),(N/L)(y + 5)}VMN/L)x ,j),((N/L)x ,(N/L) yo )) 



/(xo, B o),«, £ W; E je{W L), ,iv } Ve(((N/L)x ,j), ((N/L)x , (N/L)y )) V > V ° 



/•cast / \ , ■ ■ _ : < ; " 1 1 ■ _ ' f nr v. ~ 

•/(* ,?/o)AeW V ...nJ(i(N/T,\-, ln \((N/T,\-r. n (N/T.\i ln \\ ^ 



E t g{(Ar/L)( a: -a),(w/-L)( a: +a)} ^((», ( N / L )yo)i ((N/L)x , (N/L)y )) 



HG{(N/L)x ,N} 



rk((i, (N/L)y ), ((N/L)x , (N/L)y )) 



We consider the limits of the above four functions for e — >■ and for S — > and get the probability 
density functions of the destination distribution conditioned w.r.t. the four segments. 

«$<»> - J5sES3«S*.<») - ^ ("I 
«t) <*> = J™ ss ^<-»>.<..« = r 2 ^ (18 » 

From the above equations, it easy to see that the stationary destination distribution is uniform over 
each of the four segments. We are now able to calculate the probability that an agent, visiting point 
(#0,2/0)5 has destination in one of the four segments. 

^(xo^.e = E K((N/L)x ,(N/L)y )((N/L)x ,j) 

je{0,(N/L)y o } 

4>^xt S yo),e = E K((N/L)x ,(N/L)y )(i,(N/L)y ) 

i£{0,(N/L)x o } 

CS,e= E K((N/L)x ,(N/L)y )((N/L)x ,j) 

je{(N/L)y ,N} 

Co S to), e = E K((N/L)x ,(N/L)y )(i,(N/L)y ) 

ie{(N/L)x Q ,L} 

By taking the limits for e — > 0, we get, respectively, 

(19) 
(20) 
(21) 
(22) 



/south 


= lim ( 


/south 


Vo{L 


- Vo) 






4L(x + j/o) 


-4(4 




/west 


= lim < 


iwest 


x (L 


- x ) 




^(^o,yo) 




4:L(x + y ) 


-4(4 


+ y 2 ) 


/north 


= lim < 


-north 


yo{L 


-yo) 








4L(x + y ) 


-A{xl 


+ y 2 ) 


-.east 


= lim < 


/east 


x (L 


-x ) 




^o.yo) 


"(xo,yo),e 


4L(x + y ) 


-A{xl 


+ y 2 ) 
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G Modular MTM: Stationary Distributions 

The next proposition provides the formulas for the stationary distributions of general modular MTM. 
Such formulas will be used to prove Proposition [7] 

Proposition 15 Letm = {B,TZ) be a Route System such that the associated MTM V[fR] = (T[9t], 

is strongly connected. Let a be the stationary distribution of the Kernel o/2?[9t]. Then, the stationary 

spatial distribution s o/Z?[*H] is, for every u £ S, 

1 J see 1 1 weP(T[JH]) 1 1 P.e7?. ro 1 1 Bee 1 1 iueP(T[£R]) 1 1 Ren w 1 1 

The stationary destination distributions d o/P[SH] are, /or every u,v G 5, 
C 1 1 \ - #B,u \ - o-(w) \ - 

uW ~ s(u)A\m] ^ \b\ 2-* ircj ^ lizi 

v ; 1 J see 1 1 «,eP(T[w]) 1 1 fl.erc„ Afl c „ d =u 1 1 

Proof. Since Z?[*H] is strongly connected, from Theorem |H Z?[*H] has e unique stationary distribution. 
Let 7T be the stationary distribution of Z?[*H] and let cr be the (unique) stationary distribution of its 
Kernel. Firstly, we prove (i). From the definitions and Theorem [1] (a) , it holds that 

«(«) = E #T,u7T«T,0)) 

TGT[!H] U 

= E # T <^ V^A 77T<r( T start)m}T etart (T) 

^ ' TeT[5H] u 



To evaluate A [51] we need the following claim 
Claim 1 ET e pl r l = l Bundle (^)IE Be P^- 
Proof. Let P = (B u B 2 . . . Bfi, then 



TeP TeBt-B 



Em = E i t i 

reB 1 B 2 -...Bi 

E E ••• E(n + n + ---N) 

s 2 Si ePi 



si£Bi s 2 £B 2 Si£Pi 

i i.in' 

i=l sEBi 1 1 

Y> |s| • Bundle(P) 

; 



|Bundle(P)|Er^7 EN 

z=l 1 j| s£B, 
- 1 

|Bundle(P)|^— #P 4 

|Bundle(P)| ^ 

Be? 1 1 
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Apt] = °H^ m {w) 

weP(T[9\]) 

E < w ) E mmuT) 

weP(T[m.]) Ter[m]°^(w) 

- E «-) E iSE^E # " 



I^J ^ I-RI ^ IBundle(P) 

weP(T[m]) Ter[SK] out (i«) 1 1 Ren w 1 1 Pefl 1 y Jl 

^ \Tl w \ ^ \R\ ^ Bundle(P) ^ 1 1 WP ' T 
weP(T[va]) 1 1 Reiz^ 1 1 Pen, 1 y n reT[!K] out («0 

^) j_ i 

^ ftj ^ R ^ Bundle(P) ^ 1 1 

weP(T[<R]) 1 1 -Re7J ro 1 1 PER 1 V n TEP 

^ ctH ^ 1 ^ |Bundle(P)| ^ #B . m 

MieP(T[M]) 1 1 Retz^ 1 1 PeP 1 v n bep 1 1 

\ - o{w) J_ \p \p #pp ■ #B 

^ \n w \ ^ \r\ ^ ^ \b\ 

w£P(TM) 1 1 RETZ m 1 1 PERBEB 1 1 



^ LB I ^ IftJ ^ |P| ^ ^ P ' B 
see 1 1 ™eP(TpR]) 1 1 PeP„ 1 1 PeP 



I-BI ^ IftJ ^ |P| 
Returning to Eq. [23 we firstly observe that the following claim holds 
Claim 2 J2tep #t, u = |Bundle(P)| £ BeP iff- 
Proof. Similar to the proof of Claim Q] 
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Thus have that 

Apt] 
1 



TeTpi] u 

^ E E #r,„^)^[9t]«(T) 

weP(T[<K\) TeT[w]° ut (u))„ 



Aim] ^ L> ,u \Tl w \ ^ \R\ ^ |Bundle(P)| 

L J weP(TpR]) TeT[3*]° ut (u>) u PeTC„ 1 1 Pes. 1 v n 

= — !— v v — v 1 v & & 

Aim] ^ \K W \ ^ \R\ ^ |Bundle(P)| ^ WT ' U WRT 

L 1 wEP(T[m.]) 1 1 PeP™ 1 1 PER 1 v n T£T[<R]° llt (w) u 

1 \ - cr(w) \ - J_ \ - 1 \ - „ 

A[9t] ^ \KJ ^ \R\ ^ |Bundle(P)| ^ 7f7> 
L J ioeP(r[s«]) 1 1 Ren m 1 1 per 1 v yi rep 

1 a{w) 1 |Bundle(P)| ^ # B ,u n n] ■ r^y 

Ap] E ro E p E |Bundle(P)| E "T^f Ov Claim© 

L J weP(T[s«]) 1 1 P.eP m 1 1 Pes 1 v n bep 1 1 

_ 1 \ - <t(w) \ - 1 \ - #P,P ■ #S,n 

Apt] ^ IftJ ^ \R\ ^ \B\ 

L J wep(T[w.]) 1 1 PeP„ 1 1 see 1 1 

1 \ - #B,u \ - O-(w) \ - #P,S 

A[9tl ^ \B\ ^ \R W \ ^ \R\ 
1 1 see 1 1 weP(T[m]) 1 1 PeP„ 1 1 

This proves formula (i). 

Now consider formula (ii) for the stationary destination distributions. From the definitions, Theo- 
rem [1] (a), and very similar calculations as those used for the stationary spatial distribution, we have 

M«) = "sfu) £ #i>-7r«T,0» 

y ' TeT[K]„ AT e „ d =u 

= sUMM] E #T,u<T(T start )il;[<R}T start (T) 

! J2 E #T,u<rMmuT) 



2^ 2^ : " it?. . i Z^ i pi Z^ 



-PT 



s(ii)A[JH] ^ ^ ' \n w \ ^ \R\ ^ IBundle(P) 

0(«)Apt] Z- IftJ Z. |PJ Z^ |Bundle(P)| Z. ffT ' u ^ P ' T 

1 x a(w) 



E ilEiE IRnnrLpM E (24) 



s(w)Af5Kl ^ 17^1 ^ IPI ^ IBundle(P) 

w 1 J t«eP(r[9t]) 1 1 PeTC lu 1 1 Pes 1 y TeP^T and =v 

Notice that if R 6 1Z W is such that R en d ^ v then 

E |Bundle(P)| E #r, u = 

PeP 1 y Jl T£P/\T end =v 

Thus, returning to Eq. [Ml we have 
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MV) ~ s(u)A[m] Z- WJ L> \R\^ |Bundle(P)| ^ #7> 
1 x a(w) 



\ - o{w) x - 1 1 >jp 

Z^ II?.... I Z^ I p| Z^ |R,,nHWPM Z^ WT ' U 



s(u)hM\ ^ \1Z W \ ^ \R\ ^ |Bundle(P)| 

K ' 1 1 weP(T[<K\) 1 1 Ren m f\R snd =v 1 1 PeR 1 y n tep 

1 v-^ CT H 1 |Bundle(P)| ^ . 

^)AP] E ^ E p E |Bundle(P)| E ^ 0* Claim© 

K ' 1 1 weP(T[m]) 1 1 Ren w AR end =v 1 1 per 1 v /! bgp 1 1 

<j(w) 1 #R,b • 



Z^ hp I Z^ i pi Z^ 



s(w)Af5Kl ^ \K W \ ^ \R\ ^ LB I 



uroti Z^ iri Z^ h?..i Z^ 



s(u)Af5Kl ^ \B\ ^ \H W \ ^ \R\ 



□ 



G.l Proof of Proposition [7] 

We need, the following preliminary lemmas. 

Lemma 16 Let T> = (T, \&) be a MTM such that a function f : P(T) — > M. exists satisfying the 
following properties: 

(i) 

E f(u) > and Vu E P(T) f(u) > 
usp(T) 

(n) _ 

VuEP(T) Yl f( T ^rt)i>T start (T) = /(«) 
Ter in (u) 

Then, the following map 

W E P(T) a(u) = - f{u) 

2ZveP(T) f( v > 

is a stationary distribution of the Kernel of T> . 

Proof. From Property (a), a is a probability distribution over P{T). It holds that, for every u E P(T), 
E a(v)Pr(v^u) = E MT) 

vGP(T) veP(T) T£T(v,u) 

= E v m T ( \ E MT) 

vep{T) ^we-pcn j y w > t 

E E fWMT) 

v£P(T) T£T(v,u) 

E /( r s *-*)V-T storf (r) 

Ter in (u) 

/(m) (by Property (b)) 



J2weP(T) 




1 






/H 


1 








a{u) 





□ 
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Lemma 17 Let *H = (B, 1Z) be a balanced Route System. Then, the map a so defined 

V« e P{TM) a{u) = & 
is a stationary distribution of the Kernel of the associated MTM Z?[*H]. 

Proof. Consider the function / : P(T[fK]) -> M such that f(u) = \K U \, for every u e P(T[5H]). Clearly, 
/ satisfies Property (a) of Lemma [TBI To verify that it also satisfies Property (b) w.r.t. the MTM 
2?[9t], consider the following 

E E E ^E^y, 

TGT[SR] in (u) T6T[K] in (M) 1 ,tortl R6Rr start 1 PGP 1 V ;l 



X! X] X] i pi X] 



# 



P,T 



\R\ ^ IBundlefP) 
X X Tr[ X |Bundle(P)| X # P - T 

vGP(T[m]) R£TZ V 1 1 PGP 1 v ; TeT[<K](„,„) 

X X Tri X |Bundle(P)| X # P - T 
vEP(T[<R]) PGTC„, U 1 1 PGP 1 v Ter[SH] ( „,„) 

- x x x i^dWi Bundie ( p )' 

ti£P(T[5i]) fleK»,„ 1 1 PGP 1 v ;l 

^ |{i?eft„|7e e „ d = u}| 

i;£P(T[9t]) 

= |{JjG%|i2 en( j=tt}| 

= (since *K is balanced) 

- /(«) 

Hence, we can apply Lemma [T6l and so we obtain that the map <r, defined as 

[ ' E.GP(r)/H W 
is a stationary distribution of the Kernel of the MTM P[$H]. □ 

Proof of Proposition \7\ Since 9\ is balanced and T>[D\] is strongly connected, from Lemma IT7l the 
unique stationary distribution of the Kernel of P[9t] is the following 



We firstly consider the spatial stationary distribution. From Proposition[T5] (i) and the above expression 
for a we obtain 

/ \ _ 1 \ - #p,n \ - ajw) \ - #p : p 

[U > ~ AM] ^ \B\ ^ YR, W \ ^ \R\ 

AMI ^ TbT ^ 1-72.1 ^ "JrT 

1 1 #B,u \ -> #r,b ro _, 

2^ TrT 2^ TrT (2b > 



L J 1 1 BGS 1 1 PGK 1 
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Moreover, it holds that 



see 1 1 ™eP(rpR]) 1 1 seTC„ 1 1 



By combining Equalities [55] and [211 we obtain 



IBI 2_^ i«i 

1 1 weP(T[m]) 1 w| 

|B| 2^ 1-72.1 2_, \r\ 

1 \ " #-5 \ "* #S,S 

\K\ 2^ |B| 2^ 

1 1 see 1 1 Ren 1 1 

= -^A b pK] (26) 



#S,S 



A b [0t] ^ IBI ^ |i?| 
DL J see 1 1 i?eK 1 1 

The proof of formula (ii) is straightforward and it is very similar to the one of formula (i) . □ 



H The DownTown Formulas 

In this section, we give the spatial distribution formulas for some further kinds of cell of the DownTown 
Model. 

- Start-bundles. We consider a starting-bundle Bg~^(i,j) into a horizontal block (i.e. even i and odd 
j). Observe that the number of bundle-paths containing this bundle equlas the number of parking cells 
located into blocks whose horizontal coordinate is not smaller than j (but the cells into block 
Then, 

<r-B£j;(i,3) = m((n+l)(n-j + l)-2) (27) 
Similarly, for a starting-bundle Bti(i,j) into a horizontal block, we get 

a-B s ~ k (i,j) = m ((n+l)(j + l)-2) (28) 

As for cr-Bg^ (i, j), it is equal to number of parking cells located into blocks whose horizontal coordinate 
is less than j. So, we get 

v-Bt; k {i,j) = m((n + l)j-l) (29) 
Similarly, for a starting-bundle Bgt(i,j) into a horizontal block, we get 

a-B-+(i,j) = m ((n + l)(n-j)-l) (30) 

As for starting-bundles into vertical blocks (i.e. odd i and even j), we can get the correct formulas for 
a-Bg* k (i, j) by replacing index j with i. 

- End-bundles. We consider an end-bundle B^\{i,j) into a horizontal block (i.e. even i and odd j). 
Observe that the number of bundle-paths containing this bundle equals the number of parking cells 
located into blocks whose horizontal coordinate is smaller than j plus the number of positive parking 
stripes located into blocks whose horizontal coordinate is j (but the positive cells into block 
Then, 

a-B++(i,j) = m((n+l)j + n/2-l) (31) 
Similarly, for a starting-bundle B~^~ k (i,j) into a horizontal block we get 
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°- B s,k(h.j) = m((n + l)(n-j)+n/2-l) 
Similarly, for the other two cases we get 

°-B^; k (i,3) = m((n + l)(n-j)+n/2-l) 



(32) 



(33) 



^-^(U) = m((n + l)i + n/2-l) (34) 

As for ending-bundles into vertical blocks (i.e. odd i and even j), we can get the correct formulas for 
a-B*g k {i,j) by replacing index j with i. 

- Transit-bundles. By applying counting arguments, we get the correct formulas for transit-bundles. 
Define function 



2 (n — j + l)(nj — (j + 1)) + m 2 (n — j — l)(nj + 2j + n — i) i even, j odd with 



Then, we get 



m 2 (n - j + l)(nj - |(j + 1)) + m 2 (n - j - l)(nj + 2j + n) 

n+l , 
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\n - j + l)(nj - sfi(; + 1)) + m 2 (n - j - l)(nj + 2j) 







{rj(i,j) i even, j odd 
T](j, i) i odd, j even 
otherwise 

By exploiting symmetric properties, from the above formula we get 



i £ {0,n} 
i = and j odd 
i = n and j odd 
otherwise 



(35) 



<r-B T (iJ) = a-B^(n- 
- Cross-bundles. By applying counting arguments, we get the correct formulas for cross-bundles. 

m 2 (n 2 — n + (2n — l)j — (n — \)i — 2ij) i and j even with 



(36) 







i ^ n and j / 
otherwise 



By exploiting symmetric properties, from the above formula we get 



(37) 



-B. 



v,+- 
c 

if,-- 
c 

lV,- + ( S 



{iJ) 



a-B"' ++ (n- 

a-B"' ++ (n - i,n-j) 
h,++, 



a-B v d - + {i,j) = a-B^ ++ (j,n-i 



It holds that 



a-B%'+(i,j) = { 



m 2 (n — j)(2nj — i + j) i and j even with 

i,j i {0,n} 

m 2 (n — j)(2nj + j/2) i = and j even with 



m 2 (n-j)(|j(n+l)-n) 







n and j even with 



j<£{0,n} 
otherwise 



(38) 
(39) 
(40) 



(41) 



By exploiting symmetric properties, from the above formula we get 
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It holds that 



<• (hi) = (T-B c ' + (n - 

(n - i,n-j) 
(j, n - i) 



a-B, 



a-B 
a-B. 



H,+ i 
C 

H,+ i 
C 
H. 

C 



a-B^' + = ( m 2 (nj + n + j) i — n and j even with 



m 2 (2ij + i + j) i and j even with 
j ^ and i g {0, n} 







3^0 
otherwise 



(42) 
(43) 
(44) 



(45) 



By exploiting symmetric properties, from the above formula we get 



a-B% + + 




= °~B% A 


(n - j, i) 


(46) 






= °-B" A 


(n — i,n — j) 


(47) 


a-B v d - 




= o-B^ 


~ U, n-i) 


(48) 



Computing the Downt-Town Spatial Distribution: Relevant locations. We now use the above 
values to compute the stationary spatial distribution given by Prop. [7] We give such formulas for some 
kind of cells (the others can be obtained by the same counting arguments). Let A = K\^P]/m 2 be 
the normalization constant with AbpK 2 '] defined in Prop [3 

- [Transit Blocks]. Let u be a transit cell of index k (with any k £ {1, . . . , m}) in the positive transit 
stripe of the horizontal block (i.e. i even and j odd), (with i ^ {0, n}). Then 

Thus, since a-Bg'* h (i, j) and a-B^, + h (i^ j) do not depend on h, we have 

*(«) = (v-B$(i,j) + k ■ v-B++(iJ) + <r-B£*(h3) + (™ - k + l)a-B++(i,j) + a-B+; k (i,j)) 

= X^M^I m ( m ( n ~ ] + l ^ n] ~ ( J + X ) + rn ( n - 3 ~ 1 )( n ^' + 2j + n - i) + k((n + l)(n - j + 1) 

TL Th 

(n + l)j - 1 + (m - k + l)((n + l)j + - - 1) + (n + l)(n - j) + - - 1) 

slow(fc) / ., ,. fe , , .. 1 , 
A V mm 



where 



rt 

a(i,j) = (n-j)(2nj +j + n- i - 1) + (n- 2)j - - +i - 2 

TL 

b(j) = n (n+l) + --2(n+l)j 
c(j) = (n + l)(n + j) + n-3 

[Parking Cells]. Let u be a parking cell of index k (with any k £ {1, . . . , m}) in the positive 
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parking stripe of the horizontal block (i.e. i even and j odd). Then 



paUSS m((n + l)(n - j + 1) - 2 + (n + l)j - 1 + (n + l)j + n/2 - 1 + 



A b [$H»] 
+(n + l)(n- j)+n/2- 1) 
pause 2n 2 + 4n - 4 



(50) 



- [Cross- Ways] . Let u be a crossing cell (+, — ) in the crossing block (i,j) (i.e. i,j even ) with 
i,j {0,n} Then 

S(U) = J^]^-Bc' + (id) + cr-B^ + -(i,j) + a-B^-(i,j) + a-B^- + (i,j)) 

m 2 ((n — j)(2nj — i + j) + 2ij + i+ j + i(2n(n — i) — j + n — i) + 



+n 2 -n+(2n- 2j - l)(n - i) - (n - 

^ (n 2 (2j + 2i + 3) - 2n(j 2 + i 2 + j + i + 1) - (t - j) 2 + 2(ij + j + i)) (51) 
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